I_D=[-200,-180,-160,-140,-120,-100,-80,-60,-40,-20,0,20,40,60,80,100,120,140,160,180,200];
I_Q=0;

#inizializzo tutte le stringhe che mi serviranno
numero=0;
prefisso1='LAMBDA_1_';
prefisso2='LAMBDA_2_';
prefisso3='LAMBDA_3_';
curvaestrattaa='CURVA_ESTRATTA_A_';
curvaestrattab='CURVA_ESTRATTA_B_';
curvaestrattac='CURVA_ESTRATTA_C_';
fft1='FFT_1_';
fft2='FFT_2_';
fft3='FFT_3_';


#inizio il ciclo sui diversi valori di Id

for d in range(0,21,1):
	numero+=1;
	n=str(numero);

	#creazione grafici flusso fase 1
	
	flusso1=prefisso1+n;
	curva1=curvaestrattaa+n;
	flux1=fft1+n;

	EvolutiveCurve2D(name=flusso1,
               evolutivePath=EvolutivePath(parameterSet=[SetParameterFixed(paramEvol=VariationParameter['I_D'],
                                                                           currentValue=I_D[d]),
                                                         SetParameterXVariable(paramEvol=VariationParameter['ANGPOS_ROTOR'],
                                                                               limitMin=1.0,
                                                                               limitMax=91.0)]),
                formula=['LAMBDA_1'])



	#creazione fft di flusso1


	result = SpectrumAnalysis(curve=Curve2d[flusso1],
                 component='LAMBDA_1',
                 cycle=FullCycle(),
                 numberHarmonics=10,
                 computeDCComponent=NonContinue(),
                 scaleDecibel=FirstHarmonic(firstHarmonicValue=100.0),
                 nameRebuiltCurve=curva1,
                 nameFourierTransform=flux1)









	
	#creazione grafici flusso fase 2

	flusso2=prefisso2+n;
	curva2=curvaestrattab+n;
	flux2=fft2+n;

	EvolutiveCurve2D(name=flusso2,
         evolutivePath=EvolutivePath(parameterSet=[SetParameterFixed(paramEvol=VariationParameter['I_D'],
                                                                     currentValue=I_D[d]),
                                                   SetParameterXVariable(paramEvol=VariationParameter['ANGPOS_ROTOR'],
                                                                         limitMin=1.0,
                                                                         limitMax=91.0)]),
           formula=['LAMBDA_2'])



	#creazione fft flusso 2

	result = SpectrumAnalysis(curve=Curve2d[flusso2],
                 component='LAMBDA_2',
                 cycle=FullCycle(),
                 numberHarmonics=10,
                 computeDCComponent=NonContinue(),
                 scaleDecibel=FirstHarmonic(firstHarmonicValue=100.0),
                 nameRebuiltCurve=curva2,
                 nameFourierTransform=flux2)


	
	#creazione grafici flusso fase 3
	
	flusso3=prefisso3+n;
	curva3=curvaestrattac+n;
	flux3=fft3+n;

	EvolutiveCurve2D(name=flusso3,
             evolutivePath=EvolutivePath(parameterSet=[SetParameterFixed(paramEvol=VariationParameter['I_D'],
                                                                         currentValue=I_D[d]),
                                                       SetParameterXVariable(paramEvol=VariationParameter['ANGPOS_ROTOR'],
                                                                             limitMin=1.0,
                                                                             limitMax=91.0)]),
              formula=['LAMBDA_3'])	



	result = SpectrumAnalysis(curve=Curve2d[flusso3],
                 component='LAMBDA_3',
                 cycle=FullCycle(),
                 numberHarmonics=10,
                 computeDCComponent=NonContinue(),
                 scaleDecibel=FirstHarmonic(firstHarmonicValue=100.0),
                 nameRebuiltCurve=curva3,
                 nameFourierTransform=flux3)

	
	#estrazione grafici flusso fase 1	

	CurveVariation2D[flux1].exportExcel(xlsFile='fft_flusso1_iq0',
            mode='add')

	#estrazione grafici flusso fase 2	

	CurveVariation2D[flux2].exportExcel(xlsFile='fft_flusso2_iq0',
            mode='add')

	#estrazione grafici flusso fase 3	

	CurveVariation2D[flux3].exportExcel(xlsFile='fft_flusso3_iq0',
            mode='add')